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Abstract 

For a class of degenerate diffusion processes of rank 2, i.e. when 
only Poisson brackets of order one are needed to span the whole space, 
we obtain a parametrix representation of Mc Kean-Singer |MS67] 
type for the density. We therefrom derive an explicit Gaussian upper 
bound and a partial lower bound that characterize the additional 
singularity induced by the degeneracy. 

This particular representation then allows to give a local limit the- 
orem with the usual convergence rate for an associated Markov chain 
approximation. The key point is that the "weak" degeneracy allows 
to exploit the techniques first introduced in Konakov and Molchanov 
|KM85] and then developed in [KMOO] that rely on Gaussian approx- 
imations. 

Resume 

Pour une classe de processus de diffusion de rang deux, i.e. lorsque 
seuls des crochets de Poisson d'ordre un permettent d'engendrer l'espa- 
ce, nous obtenons une representation parametrix de type Mc Mean- 
Singer [MS67] de la densite. Nous en derivons une borne superieure 
Gaussienne explicite et une borne inferieure partielle qui caracterisent 
la singularity additionnelle induite par la degenerescence. 

Nous donnons ensuite un theoreme limite local pour une approx- 
imation par chaine de Markov associee. Le point crucial est que la 
faible degenerescence permet d'exploiter les techniques initialement 
introduites par Konakov and Molchanov |KM85| puis developpees 
dans |KM00] et qui reposent sur des approximations Gaussiennes. 
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1.1. Global overview. Let us consider in R , d > 1 the Markov diffusion 
process with generator 

L= \ H aij(x)dl iX .+ bi(x)d Xi . 
;,je[M] 2 ie[i,d| 

If the coefficients of L are smooth enough, say C 1 (]R d ), bounded, and the 
diffusion matrix A(x) = (a{j(x)) is uniformly elliptic (VA £ R d , (AX, A) € 
[5, for an appropriate 5 > 0) then the associated process (A^) 4 >o has a 
transition density p(t, x, y) which is the fundamental solution of the parabolic 
problem <9fp(.) = L x p(.), p(0,x,y) = 5 y (x). Of course, one also has dtp{.) = 
L* yP {.), p(0,x,y) = 5 x {y). 

Moreover, this density satisfies uniformly in t G]0, T] the following Gaus- 
sian bounds 

M_1 ( n, \ x -V\* \ <r U W M I \ X ~V\ 2 \ 

-par ex p ( v - M ^— J ^ 'ft ^ ^ ex p |^-^r J ' 

where the constant M depends on T, d, the ellipticity constant and the norms 
of the coefficients in C l (M. d ), see e.g. Aronson |Aro67| or Stroock |Str88j . 

The above estimations express the following physically obvious fact: if the 
process starts from xq € M. d , then for small t > 0, in the neighborhood of 
xo it is "almost Gaussian" with the "frozen" diffusion tensor A(xq) and the 
drift b(xo). 

The justification of this fact requires the solution of the perturbative inte- 
gral equation for p(-) (so-called Parametrix equation), where the leading term 
of the perturbation theory for p(-) is exactly the Gaussian kernel po(-) corre- 
sponding to the "frozen" coefficients at xq. For details concerning Parametrix 
equations we refer the reader to Mc Kean and Singer |MS67j . Friedman 
[Fri64] or |KM85| . 

If the matrix A(x) degenerates, but the coefficients a, b are still smooth, 
the diffusion process (Xt)t>o with generator L exists (one can use the Ito 
calculus for the direct construction of the trajectories), but has generally 
speaking no density. 

k 

Consider now generators of the form L = ^^4^ + Aq, k < d, where 

i=i 

(Ai) i£ p^ k j are first order operators (vector fields) on W 1 (or more gener- 
ally on smooth manifolds) with C°° coefficients. Sufficient conditions for the 
existence of the density can be formulated in terms of the structure of the Lie 
algebra of the vector fields on M. d , with usual linear operations and the Pois- 
son bracketing [.,.]. Namely, if the vector fields Ai,...,A k , [A, An](z,m)e[o,fc] 2 > 
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[Ai, [Am, ^4n]](z,m,n)e[o,fc] 3 J ••• span M. d then the density exists. This result is 
due to Hormander |Hor67| . see also Norris |Nor86j for a Malliavin calculus 
based probabilistic proof. Operators having the previous property are said 
to be hypoelliptic. Also, in |Hor67j . Hormander stressed that the seed of the 
idea of hypoellipticity goes back to Kolmogorov's note |Kol34| . 

A. Kolmogorov made the following important observation. Let d = 2. 
For the generator L = \d xx + axd y , a ^ 0, the solution of the associ- 
ated SDE writes (X t ,Y t ) = (x + W t ,y + a(x t + Jq W s ds)), where W is 
a standard one dimensional Brownian motion. Thus (X t ,Y t ) has two di- 
mensional Gaussian distribution with mean (xo,yo + ax^t) and covariance 

( t s£ \ 

matrix C = I at 2 a 2 t 3 ■ Note that the transition density for small t 
V "2" ~S~ J 

has higher singularity than the usual heat kernel. In Hormander's form 
L = \A\ + Aq, A\ = 8 x ,Aq = axdy so that L4i,^4o] = a dy and thus, 
A\, [Ai,Aq] have together rank 2. 

In this paper, using a parametrix approach derived from the work of McK- 
ean and Singer |MS67j . we are able to derive a Gaussian upper bound, and a 
"partial" lower bound with the two previous time scales, and an associated 
local limit theorem in the following case. 

1.2. Statement of the problem. We consider M. d x Revalued diffusion pro- 
cesses that follow the dynamics 

[ X t = x+ f t b(X s ,Y s )ds+ f a{X s ,Y s )dW s , 
(1.1) h J ° 

^Y t = y + ^ X s ds, 

where (Wt)t>o is a standard d-dimensional Brownian motion defined on some 
filtered probability space (SI, JF, (^ r t ) t > , IP) satisfying the usual assumptions. 
We assume that a is uniformly elliptic and that b,a are C 1 , bounded, Lip- 
schitz continuous functions so that there exists a unique strong solution to 

m- 

Such kind of processes appear in various applicative fields. For instance 
in mathematical finance, when dealing with Asian options, X represents the 
dynamics of the underlying asset and its integral Y is involved in the option 
Pay-off. Typically, the price of such options write K x [iP(Xt,T Yx)], where 
for the put (resp. call) option the function ip(x, y) = (x — y) + (resp. (y — x) + ), 
see [BPV01] and [TemOl] . 

The cross dependence of the dynamics of X in Y is also important when 
handling kinematic models or Hamiltonian systems. For a given Hamilton 
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function of the form H(x, y) = V(y) + ^-, where V is a potential and ^j- the 
kinetic energy of a particle with unit mass, the associated stochastic Hamil- 
tonian system would correspond to b(X s ,Y s ) = —(d y V(Y s ) + F(X S ,Y S )X S ) 
in fll.ljl . where F is a friction term. When F > natural questions arise con- 
cerning the asymptotic behavior of (X t ,Y t ), for instance the geometric con- 
vergence to equilibrium for the Langevin equation is discussed in Mattingly 
and Stuart |MS04j . numerical approximations of the invariant measures in 
Talay [Talf)2j . the case of high degree potential V is investigated in Herau 
and Nier |HN04| . Under the previous boundedness assumption on b, equa- 
tion (11,11) describes frictionless Hamiltonian systems with "almost linear" 
potential. 

Importantly, the two time-scales coming from Kolmogorov's example, and 
that we obtain for the density associated to (jl.ip . can be exploited to in- 
vestigate small time asymptotics of the previous models. For instance, for 
the Asian option, a normalization is required in the pay-off to make both 
quantities scale-homogeneous. 

As mentioned above, equation (jl.ip provides one of the simplest forms 
of degenerated processes and the previous assumptions guarantee that Hor- 
mander's theorem is satisfied taking only the first Poisson brackets between 
the vector fields associated to the drift and the diffusive part in (jl.ip . In 
a more general hypoelliptic setting, let us mention the work of Cattiaux 
|Cat90l l"Cat91] whose assumptions include the case (jl.lj) . but who obtains 
less explicit controls, see his Proposition (1.12). Under the "strong" Horman- 
der condition that involves the Poisson brackets of the diffusive part of the 
process, small time asymptotics of the density are discussed in Ben Arous 
|Ben88| or Ben Arous and Leandre |BL91j . Eventually, in whole generality 
two-sided bounds for the density of degenerate diffusions are investigated in 
Kusuoka and Stroock |KS87| . All these work strongly rely on Malliavin cal- 
culus techniques. We want to stress that the parametrix approach is not very 
well suited to study general degenerate processes. Anyhow, the counterpart 
is that it gives by construction more explicit controls. In the non-degenerate 
case, for a-H61der continuous coefficients, it directly gives two-sided Gaus- 
sian estimates. The lower bound on the diagonal in small time derives from 
the series representation and the global lower bound is obtained thanks to a 
chaining argument as in |KS87j . Here, we still derive a lower bound in small 
time from the series, but do not succeed to do a chaining argument 

Also, our controls remain valid if the coefficients in (jl.ip are uniformly 
a-H61der continuous, a case for which Hormander's Theorem breaks down, 
see Section l3l Remark 13. II for details. 

A natural question then concerns the Markov chain approximation of 
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(jl.ip . For non degenerated processes this aspect has been widely studied, 
see e.g. [KMOOj for local limit theorems. In |BT96j . using Malliavin calculus 
techniques, Bally and Talay obtain an expansion at order one w.r.t. the time 
step for the difference of the densities of the diffusion and a perturbed Euler 
scheme, i.e. the stochastic integrals are approximated by Gaussian variables 
and an artificial viscosity is added to ensure the discrete scheme has a den- 
sity. This rate corresponds to the usual "weak error" bound. Since we follow 
the local limit theorem approach we can handle a wider class of random 
variables in the approximation but obtain a rate of order 1/2 w.r.t the time 
step. Of course, plugging Gaussian random variables in our approximation 
yields to rate h as in [BT96J. 

Importantly, as opposed to |BT96| . we do not need to introduce an ar- 
tificial viscosity to ensure the existence of the density for the underlying 
degenerate Markov chain. We develop analogously to the continuous case a 
parametrix approach to express the density of the Markov chain in term 
of the density of an auxiliary frozen random walk. The random walk is 
degenerated as well, but has a density after a sufficient number of time 
steps, see Subsection 14.41 for details. The local limit theorem is then derived 
from an accurate comparison of the parametrix expansions of the densi- 
ties of the process and the chain. To motivate this result we can consider 
the case of the approximation of a "digital Asian call" i.e. of the quantity 
P[(T _1 y T - X T )+ > K] for a given K G M + . Indeed, the local limit theorem 
associated to our scheme directly relates the densities of the discrete and 
continuous objects which is not the case if we only consider a discretization 
of the non degenerate component and a numerical estimation of the integral, 
since in that case the approximating couple can fail to have a density. 

The paper is organized as follows. In Section fl~3l we give our assumptions 
and fix some notations. Then, since the form of the Markov chain approxi- 
mation strongly relies on the proof of our results for the diffusion we choose 
to divide this paper into two parts. Sections [2] and [3] deal with the results 
for the diffusion and their proofs. Section [4] is dedicated to the Markov chain 
approximation of (jl.lj) . the associated convergence results and the key points 
of the proofs. The complete proof of the local limit Theorem can be found 
in the Appendix. 

1.3. Assumptions and Notations. We suppose that the coefficients of 
equation (jl.lj) satisfy the following assumptions. 

(UE) 3(A min ,A max ) G (0,oo) 2 , Vz G R d , V(x,y) G R 2d , A min |z| 2 < 
{a(x,y)z, z) < A max |z| 2 , denoting a(x,y) = aa*(x,y). From now on we sup- 
pose that a is the unique symmetric matrix s.t. aa = a. We are interested 
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in the density of the process and its approximation at a given time. Hence, 
from the uniqueness in law, the previous assumption can be made without 
loss of generality. 

(B) The coefficients b,a in (jl.ljl are C , uniformly Lipschitz continuous 
and bounded. 

Throughout the paper we consider the running diffusion (|l.ip up to a 
fixed final time T > 0. We denote by C a generic positive constant that 
may change from line to line and only depends on T, and the parameters 
appearing in (UE), (B). We reserve the notation c for constants that only 
depend on parameters from (UE), (B). Other possible dependencies are 
explicitly indicated. 

2. Explicit parametrix and associated controls for the density 
of the diffusion. The previous assumptions guarantee that Hormander's 
Theorem, see e.g. Nualart |Nua98| . holds true, and therefore that Vi > 
0, (X t ,Y t ) has a density w.r.t. the Lebesgue measure. Introduce the vec- 
tor fields in 



»2<2 



A (x,y) 



( h(x,y) \ 
bd{x,y) 

V x d j 



Vj G [l.dj, Aj(x,y) 



( vij{x,y) \ 

°dj(x,y) 




(2.1) 



One directly derives the following 
Proposition 2.1 For all (x,y) G R 2d , 

Span(,4i(x,y), A d (x,y), [A (x, y), A^x, y)], 



[A (x,y),A d (x,y)]) 



where V(z,j) G 
bracket. 



<i] 2 , [Aj,^4j] = AiVAj — AjVAi denotes the Poisson 



Fix T > and < t < T, (x,y) G R 2d . Since, we now know that 
(X t ,Y t ) has a transition density, i.e. F[X t G dx',Y t G dy'\Xo = x,Yq = 
y] = p(t, (x, y), (x' , y'))dx'dy', our aim is to develop a parametrix for (II. ip to 
obtain an explicit representation of this density. 

Recall that we consider the following SDE 



(2.2) 




: b(X s , Y 8 )dt + a(X s , Y s )dW s , X 
X s ds, Y = y. 



x, 
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For the parametrix development we need to introduce a "frozen" diffusion 
process, (X s ,Y s ) se[0 $ below. Namely for all {x',y') G R 2d , define 




dX y,y' = a (x', y' - x'(t - s))dW s + b(x', y')ds, X 1 / > y ' = x, 

d?y>y' = xy>y'ds, y 1 /^' = y. 



The key point is that the above process is gaussian. The arguments in the 
second variable of the diffusion coefficient can seem awkward at first sight, 
it includes the transport of the frozen point x' with a time reversal. This 
particular choice is actually imposed by the natural metric of the frozen 
process, see Proposition [3TTJ in order to allow the comparison of the singular 
parts of the generators. 

The processes (X S ,Y S ) and (Xi> x '>y',Yy'y'), s G [0,t], have the following 
generators: V(x,y) G M. 2d , tp G C 2 (M 2d ), 

Lip(x,y) = QlV (a(x,y)Dlipj + (b(x,y), V x ip) + (x,V y ip)j(x,y), 
L</^(x, y) = QTr (a (x' , y> - x'(t - s)) D 2 ^) + (b (x' , y') V^> 

(2.4) +(x,V^))(x,y). 

From these operators we define for < t < T, ((x, y), (x',y')) G (R 2d ) 2 the 
kernel H by 

H(t, (x, y), (x, y')) = (L - L)p(t, (x, y), (x',y')), 

where p(t, (x, y), (•, •)) := (x,y), (•, •)), £ := ^ ,y respectively 

stand for the density of the process (xl' x ' v , Y^' x ' y ) and the generator of 
(X^ ,x < y ,Yg' x ' y ) sg [o,t] at time 0. We omit to explicitly emphasize the depen- 
dence in t, x', y' for notational convenience. 

Remark 2.1 Note carefully that in the above kernel, because of the linear 
structure of the model the most singular terms, i. e. those involving derivatives 
w.r.t. y, vanish. 

The next proposition gives the expression of the density p in terms of an 
infinite sum involving iterated convolutions of the density p with the kernel 
H. Namely, 

Proposition 2.2 (Parametrix expansion for (12. 2h ) 

For allO<t<T, ((x, y), (x', y')) G (R 2d ) 2 , 

+oo 

(2.5) p(t, (x, y), (x',y')) = £p® (t, (x, y), (x', y>)), 

r=0 
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where 

f ®g(t,(x,y),(x',y')) = [ du [ f(u,(x,y),(z,v)) 

JO JR 2d 

xg(t - u, (z,v), (x',y'))dzdv, 

= p and = H ® H^ r 1 \ r > denotes the r-fold convolution 
of the kernel H . 

The previous Proposition is a direct consequence of the usual parametrix 
recurrence relations. For the sake of completeness we provide its proof in 
Section [3l 

Now, since for < t < T (X s ,Y s ) se \ 0tt ], is a Gaussian process, p and its 
derivatives are well controlled. The previous expression is the starting point 
to derive the following 

Theorem 2.1 (Gaussian bounds) 

There exist constants c, C > s.t. for all < t < T, ((x, y), (x' , y')) £ 
(R 2d ) 2 , one has: 

p(t, (x,y), (x',y')) < Cp c {t, (x,y), (x',y')) 

(2.6) 
where 

p c (t, (x,y),(x',y')) := ^ 2 ^ xexp 

enjoys the semigroup property, i.e. VO < s < t < T, 
/ dwdzp c (s, (x,y), (w,z))p c (t - s, (w,z), {x ,y')) = p c (t, (x,y), (x',y')). 

JR 2d 

Also, for a given Co > 0, 3to := to(Co,c, C) s.t. for t < to, [ ^ ^ — h 
3 \y'-y-^\^ < p(tj (Xj y)) ^ ^ > c-%-i(t, (x, y), (a/, y'))- 

Remark 2.2 T/te /ower bound, obtained in small time and compact sets, 
derives from the parametrix representation of Proposition \2.S\ and the upper 
Gaussian control. It remains an open problem to find a well suited chaining 
argument to derive a global lower bound for this degenerate case. 

Remark 2.3 Note that the above result would remain valid if we replaced 
the dynamics of Y t in (TTT]) by Y t = y + f Q t F(X s )ds for a C 2+a , a > 0, 
Lipschitz continuous mapping F : R d -> R d s.t. VFVF* is non degenerated, 
i.e. 3c ,V(£,x) eR d xR d ,\(VFVF*(x)C,0\ > c o|£| 2 - Indeed, in such a case, 
(X s , F s ) s6 [o,t] := (F(X s ),Y s ) s£ [o,T] follows a dynamics of type ([□]). 



At 



+ 3 



y -y 



(x+x')t 



t 3 
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3. Proof of the main results: diffusion process. 



3.1. Proof of Proposition \2.Bc parametrix expansion. Following Cattiaux 
|Cat90j and Lemma 13.11 one derives that p, p have continous densities with 
bounded derivatives. Hence, from the forward and backward Kolmogorov 
equations associated to (X, Y), (X, Y) and denoting by L* the adjoint of L, 
we have 

P(t, (x, y), (x',y')) - p(t, (x, y), (x , y')) 
t o r 

du ^T / dw dzp(u,(x,y),(w,z))p(t-u,(w,z),(x',y')) 



o 



du 



du I dwdz 



dp(u, (x,y), (w,z)) 



du 



p(t - u, (w,z), (x',y')) 



+p(u, {x,y), {w,z)) x 



dp(t - u, (w,z), (x',y')) 
du 



<1u 



du 



dwdz [L*p(u, (x, y), (w, z))p(t - u, (w, z), (x',y')) 
-Lp(t - u, (w,z), (x',y'))p(u, (x,y), {w, z)) 

t dwdzp(u, (x,y), (w,z))(L - L)p(t - u, (w,z), (x',y')) 
= p <g> H(t, (x,y), (x 1 ,y')). 



A simple iteration completes the proof. 
3.2. Proof of TheoremlEJl 



□ 



3.2.1. Proof of the upper bound. The proof is divided into two parts. First 
an elementary control on the density of (X, Y) is stated in Lemma l3~Tl Then, 
this control is used to control the kernel H and the convolutions. 

Step 1: Gaussian controls for (X,Y). 

Lemma 3.1 There exist constants c > 0, C > 0, s.t. for all multi-index a, 
\a\<3,\/0<u<t<T, \/(w, z), (x',y') G R 2d 

c d 2> d / 2 



\£%p(t - u, (w, z), (x', y'))\ < C(t - u)-H/ 2 7n 



x exp — c 



(2vrt 2 ) 

2 



\x' — w\ 2 \y' — z — \{x' + w){t — u)\ 



4(t - u) 



(t - uf 



:= C(t - u)~~p c {t - u, (w, z), (x , y )), 
p(t - u, (w,z), {x ,y)) > 2C~ 1 p c -i{t - u, (w,z), (x',y')), 
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where p c is as in Theorem \2.1\ and enjoys the semi-group property. 

The proof is postponed to Section 13.2.21 

Step 2: Control of the kernel. Recall that under (B), the coefficients 
a, b are uniformly Lipschitz continuous. Hence, it is easy to get from Lemma 
13.11 and the previous definition of H that, up to a modification of c > in 
Pc that 3Ci > 0,Wi G [0,t), 

Q 

(3.1) \H(t - u, (w, z), (x',y'))\ < p c (t - u, (w, z), (x , y')). 

y/t — U 

Lemma l3~T1 also yields that 3C2 > 0, Vu G (0,t], p(u, (x, y), (w, z)) < 
C2Pc(u, (x, y), (w, z)). Setting C := C\ V C2, we finally obtain 

\p®H(t, (x,y), (x',y'))\ 

< [ du [ p(u,(x,y),(w,z))\H(t - u,(w,z),(x' ,y'))\dwdz, 
JO JR 2d 

< [ du f C 2 p c (u, (x,y), (w,z)) 2. — p c (t - u, {w,z), (x',y'))dwdz 

JO JR™ y/t — U 

<CV/ 2 i?(l,i)p c (t,(x,y),(x',y')), 

using the semigroup property of p c in the last inequality and where B(m, n) = 
Jq 1 <iuu m_1 (l — u) n ~ l denotes the /3-function. By induction in r, 

r® ffW( tj (x> y) , < CT+WBQ., \)B(\, I) x ... x B(^±l, I) 

(3.2) xp c (t,(x,y),(x' ,y')), r G N*. 

This implies that the series representing the density p(t. (x, y). (x'. y')) 

00 

p(t, (X, y), (x', y')) = £ 2/), (*',*/')) 

r=0 

is absolutely convergent and the following estimate holds 

\p(t,(x,y),(x',y'))\ < Cp c (t,(x,y),(x',y')). 



□ 



Remark 3.1 Note carefully that the above series still converges if the coef- 
ficients b, a are only uniformly a-Holder continuous. In such case Horman- 
der's theorem does not hold, but one can show by standard techniques, see 



e.g. Baldi jBd7^ , that p(t, (x,y), (., .)) := YjP <8> H {r > (i, (x, y), (., .)) 



is a 



probability density and derive with a Dynkin like argument, see e.g. Theo- 
rem 2.3 in \Dyn63y , that it corresponds to the density of the weak solution of 
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3.2.2. Proof of the partial lower bound. From the previous proof and the 
gaussian nature of (X t ,Y t ), see Lemma I37T1 one gets 

P(t, (x,y), {x ,y')) > p(t, (x,y), (x',y')) - Ct 1/2 p c (t, (x,y), {x ,y')) 

> 2C- 1 p c - 1 (t, (x, y), (x',y')) - Ct l ' 2 p c {t, (x, y), (a/, y 1 )) 

> C~ x p c -i{t, (x,y), (x',y')) 

for + 3 \y'-y-^'+*W 2 < Cq and t small enough . 



3.2.3. Proof of the technical Lemmas. 

Proof of Lemma 13.11 We prove the statement for \a\ = 0, i.e. without 
derivation. Indeed, since our computations only involve a finite number of 
derivations that introduce some polynomials in front of the exponential, the 
general bound can be derived similarly and the result holds taking the worst 
constants. Also, with respect to the statement of the lemma, we suppose 
w.l.o.g. u = for notational convenience. We get from (|2.3p with x = w, y = 
z that for all < t < T, 

% = ^z + wt + b(x',y')^\ + J* y\{x',y' -x'(t-s))dW s }dv 
:= 1712+ + A t , 



ft 

At 



f (t - s)a (x',y' - x'{t - s)) dW s := f (t - s)a s dW s 
o J o 



(3-3) 



using Ito's formula for the last equality. Setting Vs G [0, t], a s = a 2 , recall 
from (UE) that a s is symmetric, we finally obtain that the covariance matrix 
Et of the vector (X t , Y t ) is equal to 

s _ ( Ioa s ds fo(t-s)a a ds \ 
1 V /o(* ~ s)a s ds fv(t - s) 2 a s ds J ' 

Note from (UE) that: 3c > 0, Vs G [0,T], V£ G R d , (o a f,f) > c|£| 2 . Hence, 
by the Cauchy Schwarz inequality 

2d /vi. t: z\ -> ninm.7. 7\ r. — ( tId 2 Id ^ 



VZ G K , (E t Z, Z) > c/2(C t Z, Z), C t : 



2 J d 3 J ( 



where Ct is actually the covariance matrix of a d- dimensional standard Brow- 
nian motion and its integral. 
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The mean vector of (X t ,Y t ) is equal to (ro^^m^t), with mi i( = w + 
b(x',y')t and m,2,t as in (13, 3p . Note that Ct = T.AT*, where 

r * = ( h \h\ A (tl d \ 



o / d I gl d 



Hence, C-^tn-U-T-^- 1" J ( „ I A J ( -ft 4 
Now, VZ e M M , £ := - (E^Z, z) < — c/2 (A(T~ 1 Z), T~ 1 Z). In particular, 
for Z = (Z 1 ,Z 2 ),Z 1 = x' - {w + b(x',y')t), Z 2 = y' - (z + tut + b(x' 

, m_i w / x' — w — b(x'y')t I „. . , . 

we get i Z = , i , , \ . We therefore derive 

B \ y' - z - + w)t J 

c ^ C I ' uf I /\j.|2 6c. , 1., . |2 

£ < F - io - 6(x ,y )t\ - -p\y - z - -(x + w)t\ . 

From (B) (boundedness of 6), we derive that there exist c, C > s.t. 

£ < C- 



w | 2 +3 \y' -z-±(w + x>)t)\ 2 



At t 3 



Eventually 

Cc d 3 d / 2 

p(t,(w,z),(x',y')) < ex P 



k'-^l 2 , Jy' -z- \{w + x')t)^ 



+ 3- 



At t 3 
:= Cp c (t, (w,.z), {x',y')). 

Note from |Kol34j that p c enjoys the semigroup property. This gives the 
statement for |a| = 0. The lower bound is derived similarly from the control 
VZ £ R 2d , (E t Z, Z) < ^-(C t Z, Z) achieved for c small enough. □ 

4. Markov Chain approximation and associated convergence re- 
sults. 

4.1. Global strategy. Let us recall the strategy to derive a local limit the- 
orem for the Markov chain approximation associated to a diffusion process. 
Suppose the underlying diffusion has a density with parametrix representa- 
tion as in Proposition 12.21 If the "natural" Markov chain associated to the 
diffusion has a density, the main idea is to introduce a Markov chain with 
frozen coefficients that also has a density so that the density of the Markov 
chain can be written in parametrix form as well with a suitable discrete 
kernel. 

The next step consists in comparing these two parametrix representations. 
To this end, two key steps are needed: 
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1. The comparison of the densities of the frozen Markov chain and frozen 
diffusion process. 

2. The comparison of the kernels. 

The first step relies on Edgeworth like expansions, see e.g. Bhattacharya and 
Rao |BR76j . the second one on careful Taylor like expansions. 

The local limit Theorem is then derived by controlling the iterated con- 
volutions of differences of the kernels. This procedure has been applied suc- 
cessfully in [KMOOj to derive a local limit theorem for the Markov chain 
approximation of uniformly elliptic diffusions with bounded coefficients. 

In our current framework new difficulties arise. First of all it is not obvious 
to derive that a "natural" Markov chain associated to (II. ip has a density. 
To guarantee such an existence a common trick in the literature consists 
in adding an artificial viscosity term in the discretization scheme, see e.g. 
|BT96j . Our strate gy is here different. Namely, we manage to obtain a den- 
sity for the natural frozen Markov chain deriving from (|2.3p after a sufficient 
number of time steps. We therefore consider a "macro scale" frozen model 
corresponding to this number of time steps. We then obtain a good com- 
parison between the densities of the "aggregated" chain at macro scale and 
the frozen diffusion process. This first step gives the structure of the random 
variables involved in the approximation in order to have the comparison of 
the densities. These variables have a density. From these variables, we then 
derive the Markov chain dynamics by letting the coefficients vary at macro 
scale. 

A second difficulty is that the second component in (II. ip is unbounded. 
This yields to handle a supplementary term w.r.t. the analysis carried out 
in |KM00| and to a slightly different version of the local limit theorem. In 
the sequel we first give the dynamics of the Markov chains at macro scale 
and state the local limit Theorem (Subsection 14. 2|) . We give the Lemma for 
the comparison of the densities (Subsection 14.31) and prove the existence of 
the density for the aggregated "frozen" Markov chain (Subsection 14. 4|) . The 
whole proof of the local limit Theorem is carried out in the appendix. 

4.2. Models and results. Now, fix T > 0, N G N* and let h = T/N be 
the "micro" time discretization step. Let n G N* be large enough so that the 
natural "frozen" chain associated to (|2.3|) has a density, see Proposition 14.21 
and define the "macro" scale time step h = nh and set N = N/n G N* the 
total number of "macro" time steps over [0, T]. 

For all i G [0, NJ set tj := ih. For any (x, y) G R 2d , we define on the time 
grid {0, ...,t N ] an R 2d valued Markov chain (Z£) i6[0jJV] = ((X%,Yg)*) ie[0iN] 
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whose dynamics is given by 

Z$ = (x,y)*, and Vz G [0, iV — 1], 

Y t h i+1 = Y>: + (X^ + ^b(Z^h + a(Z t h jVh^ +1 )h, 

(4.1) 

where 7„ := (1 + ^). The variables (tfi);e(o,7V] := (rjj, r]f) i( z^,Nj are i.i.d. cen- 
tered 2d-dimensional random variables. The density </n(?7i)?72) of i?i satisfies 

(Al) Eft] = 0, and Cov(^) = ( x ^ff d ) . 

(A2) There exist a positive integer S' and a function tp : M? d — > R with 
sup MgR 2d < 00 an d / ||w|| ip{u)du < oo for S = AdS' + 4 such that 

\D u u q n {u)\ <t/,(u) 

for all \u\ G [0,4]. The main result of the section, i.e. Theorem 14.11 is stated 
in terms of S' . 

We finally need a "frozen" time inhomogeneous Markov chain. For (x, y), (x', y') G 
R 2d , j G (0,iV] 2 we define (Z£) ie[0)j ] = ((X£,Y$)*) ieloJi by 

(x,y)*, and Vi G [0,j — 1], 
The i.i.d. variables (fjj , fff )ie(o,jj have density q n {.). 

Remark 4.1 Note that the models introduced in (|4.ip and (14. 2p can seem 
awkward at first sight. They actually derive from computations that yield the 
existence of the density for the natural frozen Markov chain associated to 
(|2.3p after n "micro" time steps h, i.e at the "macro" level with time step h. 
This is developed in Subsection \4.4\ 

From now on, p h {tf, (x,y), (•,•)) andp^' ' (tf, (x, y), (•, •)) ■= Ph(tj',(x,y), 
(•,•)) denote the transition densities between times and tj> < tj of the 
Markov chain (14. 1|) and "frozen" Markov chain (|4.2I) respectively. Introducing 
a discrete "analogue" to the generators we derive from the Markov property 
a relation similar to (|2.5p between ph and ph- 



Z, 



X 



Y, 



o 

h 

ti+l 

h 

U+i 



(4.2) 
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For a sufficiently smooth function /, define and by 

Lhf(tj, (x,y), (x',y')) = 



h- 1 



Ph(h, (x,y), (u,v))f(tj - h, (u,v), (x',y'))dudv 
- f(tj - h, (x,y), (x',y')) 



h- 1 



L h f(tj, (x,y), (x',y')) = 
Ph' X ' V {h, (x,y), (u,v))f(tj - h, (u,v), (x',y'))dudv 

- f(tj ~ h, (x,y), (x',y')) 



Note that because of technical reasons, there is a shift in time in the above 
definitions, i.e. the time is tj — h, instead of the "expected" tj, in the right 
hand side of the previous equations. 

A discrete analogue Hh of the kernel H is defined as 

H h (tj,(x,y),(x',y')) = (L h -L h )p h (tj,(x,y),(x',y')), < j < N. 

From the previous definition 

H h (tj,{x,y),(x',y')) = hT 1 x 
Ph -p h J ' x ' v (h, (x,y), (u,v))p h 3,x ' y (tj - h, (u,v), (x ,y'))dudv. 

Analogously to Lemma 3.6 in |KM00| we obtain the following result. 

Proposition 4.1 (Parametrix for Markov chain) . 

Assume (UE), (B) are in force. Then, for < tj < T, 

j 

(4.3) p h (tj, (x, y), (x',y')) = ]T (ph ®h H^) (tj, (x, y), (x',y')), 

r=0 

where the discrete time convolution type operator <S>h * s defined by 

(g ®h f)(tj, (x,y), (x' ,y')) 

j-i . 

= ^2 h / 9(ti,(x,y),(u,v))f(tj -ti,(u,v),(x',y'))dudv, 
i=o J 

~ (0) (r) (r— 1) 

Ph <S>h H h = ph and H h = Hh <S)h H h denotes the r-fold discrete con- 
volution of the kernel Hh- W.r.t. the above definition, we use the convention 
thatp h ® h H { h r \0, (x,y),(x',y')) = 0,r> 1. 
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Now (14. 3ft and (12. 5p have the same form. Comparing these two expressions 
we obtain the following local limit Theorem. 

Theorem 4.1 (Local limit Theorem for the densities) . 

Assume (UE), (B), (A-l), (A-2) hold true. Then, 3c > 0, 



sup 

(x,y),(x',y')£R 2d 



(l + \x'\ + \x\) sup p c (T{l + 5),{x,y),{x',y')) 

56 [0,1] 

/ / / , \ \ -i -1 

( , i r fx +x 

+xvt \ x ~ x >y -y~ T 



2 

x \p h (T, (x, y), (x',y')) - p(T, (x, y), (x',y'))\ = 0{h 1 ' 2 ), 

where p c is as in Theorem \2.1\ ph denotes the density of the Markov chain 
KTh and V(p, u, v) 6 R + x R 2d , 

X P (u,v) = p- M X (u/p,v/p 3 ), x (u,v) = (l + dnp + H 2 ) 5 '- 1 )" 1 . 

Note from the above result that the bigger is S", the better is the control on 
the tails. 

4.3. Comparison of the discrete and continuous frozen densities. The 
first step for the error analysis is achieved with the following 

Lemma 4.1 There exists C > 0, s.t. for all j G (0, Nj, p 2 := tj, 

\(Ph-p)(tj, (x,y), (x',y'))\ < Ch^-^p^Cpix' -x,y' -y- X ^ - p 2 ), 
(4.4) 

where C p {u,v) = p~ Ad ((u/ p,v/ p 3 ),((u,v) = 1+[|u|2+| * |2]( s-4)/2 , S hem 9 in- 
troduced in (A2). 

Proof. Iterating (14. 2ft from till tj we get 

1 j ~ l 

X\ = x + b{x',y')p 2 + p {—J2<T(x',y' -x'(p 2 -t k ))rjl +1 } 

3 k=0 

Y t h = y + X p 2 + £b(xi,y>)(l + ±) 
3 2 nj 



I 1 i ~ 1 1 

3 k=0 3 



(4.5) + 4j2 ZX 1 '- y - x V - f 1 - — 
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Introduce 



rrii 



and 



e 



3 ■- 



x + b(x', y')p 2 
y + xp 2 + ^-6(x',y')7n, 



??? 



m 



ln,j ■= 1 + 



3-1 



fc=0 
3-1 



1 



fc=0 



3-1 



+ 7172 E^'- - x V - ^M+i 1 



fc=0 



J 

k + 1 

3 



The dynamics of (14. 2|) thus writes 



nj 



X 1 

3 



pldxd Odxd 
Uxd 



Odxd /0 3 I a 



©3- 



Setting Vs € [0,p ], 0(s) := inf{tj := ih : t{ < s < ti + i},a s := a (x',y 
x'(p 2 - s)) we get Vj := Cov(@j) = 



1 f*. 
77 



h So dsa <t>(s) F J? dsa ^) F i h {^( s )) 

3 

J Q h dsa^ a )Ff' h (4>{8)) -p Jo dsa Hs) Fi th ((f)(s)) 



where ffW)) := + - (</>(*) + fc)M' fc (0(*)) := + &) + 

7 „/i(*j - + + (t; - (<l>(8) + h)) 2 ]. 

Thus, for h small enough, the covariance matrix Vj is uniformly invertible 
w.r.t. the parameters n,j, € N*. Denoting by g n the density of the normalized 



— 1/2 

sum 

Vj 9 3 

we derive 



p h (tj,(x,y), (x',y')) 



p 4d det(V- /2 ) 



9n V 



-1/2 



x' — m 1 - 
1 

, P 2 
y — "J^ 



Applying the Edgeworth expansion for g n (the key tool is the normal approx- 
imation of Bhattacharya and Rao, Theorem 19.3 in |BR76j ) and exploiting 
arguments similar to those of the proof of Lemma 13.11 we obtain 



Ph(tj, (x,y), (x',y')) 



1 



■go ^ 



-1/2 



x'—m 1 - 
j_ 

, P 2 
y'—vnq 



(4.6) 



p 4rf det(y/ /2 ) 
< ChV 2 p-\ p {x< - x, y> - y - ^p 2 ), 



18 



V. KONAKOV ET AL. 



where go stands for the standard 2d dimensional Gaussian density. To con- 
clude the proof, recall from the proof of Lemma 13.11 that 

P(tj, (x,y), (x',y')) 



p 4d det(Cj /2 
(4.7) 




x + b(x',y')p 2 \ ( ml 
where m c ,j = 2 4 : = 2 , and G 

' y + xp 2 + ^b(x',y') I \ m 2 c ,j 



c.j 



j- Jq 3 dsa s i $ dsa s (tj - s) 

^ 3 

4 J J dsa s (tj - s) 4 / J dsa s (tj - s) ; 



The result eventually follows from (14. 6p . (14. 7 j) and standard computations 
involving the mean value theorem. □ 

4.4. Existence of the density for the aggregated frozen process. Let ho > 
be a given fixed time step. For i £ N set U := iho- Fix (x',y') £ R 2d ,t > 0. 
We consider the frozen model defined by Xq° = x, Yq° = y and for all i E N, 

= X^ + b(x',y')h + o-(x',y' - tx')y/h^l +1 , 
Y t h ° = Y t h °+X^h 

(4.8) = Y*« + h X?° + h 2 b(x', y') + h 3 Q /2 a(x', y' - tx% +1 , 

where (£i)igN* are i.i.d., centered with identity covariance. The aim of this 
section is to show that for i large enough (X? ,Y t ho ) admits a density. We 
refer the reader to the work of Yurinski |Yur72j or Molchanov and Varchenko 

|MV77| for related topics. 

/ j^/io _ x * \ 

Conditionally to + and iterating the frozen model we get 

V Y u =y ) 

= x * + (nh )b(x',y') + a(x',y' - x't)V^h$l 

(4.9) Y t )l n =y* + (nh )x* + ^(nh ) 2 b(x', y') + (nh ) 3 / 2 a(x' ,y' - x't)$% 

where we recall 7n = (1 + ^) and 

~(i) 1 /~ ~ ~ \ 

ii : n = 7^ + £i+2 + ••• + ii+n J , 

eg = l(& fl + (i-l)& f8 + ... + (i-^)(J. 

vn V n n 



PARAMETRIX FOR SOME DEGENERATE PROCESSES 



19 



We have 



[l-^) 2 + ... + l 2 



n 



2n 2 + 3n + 1 _ 1 1 
fe2 " 3 7n(1 + 2V 



Hence, the covariance matrix of the 2d dimensional vector [QniQn) m 
non-degenerate for n > 2. 

Q n > £i n 

M 2d we derive the following 

Proposition 4.2 Le£ c/>(r) := E exp u(£i,t)J , t G M. d denote the char- 
acteristic function of the (£i)ieN* . If for all multi index (3, \/3\ = S + 2d + 1, 
\D^4>(t)\ < C(l + |r| 4+M+1 ) _1 , then for n large enough and for all multi 
index a, \a\ < 4, one has 

[ Kru^lD^+^niTuT^dnd^ < 00. 

JR d X~R d 

In particular, by Fourier inversion the density 



(4.io) ue u e 2 ) 



i 



exp(-i{(9 1 ,6 2 )*, (Tl,T 2 )*))tp n (Tl,T 2 )dTldT2 



(27r) 2rf 

exists and there exists C s.t. for all multi index u, \v\ < 4, 

C 



\D u f n (e 1 ,e 2 )\ < 



1 + 



l ( 1; 2) |s + 2d + i --MOM. 



Proof. Write 



n-l 



E 



e*p{^4i!)+<(^)}] = II 



Vi + (l-i)r 2 ' 



3=0 



<Pn(Tl,T 2 ) 
(4.11) 

We partition the space M 2d into the following disjoint sets 
A := {(r 1 ,r 2 )GlR M :|r 1 |>(l-i)|T 2 ||, 

At := ((ri,r 2 )eR 2<i :(l-^)|T 2 |<|ri|<(l--)h 
I n n 



A,-! : = ( 



;,2</ 



1 



Ji,t 2 ) G Jk : ti < - r 2 
n 
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If fa, r 2 ) G A then for % G [2, re - 2] 



ri + (1 - ^)r 2 



> 



1 



ri|-(l--)|r 2 
n 



> _ (1--)| T2 | -(!--) | T2 



n 



n 



i - 1 



\ T 2\ 



and similarly 



n+(l-i)r 2 



v 7 " 



> -^-i I-7-1 1 . Hence, 



(4.12) 



ri + (1 - ^)r 2 



2(2+1 



If (ri,r 2 ) G Af* for some i*, i* G [l,n - 2] and I G [2, re — 1 — i*] then 
elementary computations yield similarly 

2d+l 

Ti -I- I I — - — — )Tn 

(4.13) 



n + (1 - ^)r 2 



and for I G - 1] 



(4.14) 



n + (1 - ^)r 2 



/I 



2<f+l 



2 n 3d+3/2 



;2d+l 

> -U™ 1(71,71)1^ 



2 n 3d+3/2 



If (ti,t 2 ) G A n _i then for i G |1, n - 1] 

2d+l 



(4.15) 



n + (i - ^)r 2 



> 



1 



2n d+l/2 



1 



i + 1 



2d+l 



l(n,r 2 )| 



2d+l 



Use now the growth assumption on eft and the inequality 1 + J2jLi Pj < 

re pj 



Ilj!=i(l + Pj) where pj > 0, to derive from (14. lip 



\<Pn(n,T 2 )\ 



n 

J=0 



< 



ri+(l-i)re 



2d+l N 



< 



ti+(1-£)t 2 



2cZ+l 



Now equations fl4.12j) . (14,13)1 . (j4.14|) . (|4.15D yield that there exists re large 
enough s.t. 

C(re) 



bn(ri,r 2 )| < 



1 + |(ti,t 2 



\2d+l 
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where C(n) — ► +00. Anyhow, for such a fixed n, one has (p n £ L 1 (R 2d ) which 

/~fl) ~(2)\* Oft 

implies the existence of the density f n of the vectors [Q^Qn ) e ^ • The 
properties concerning the growth and derivatives of f n are derived from (|4.10l) 
and the growth and smoothness properties of <j>. □ 

Hence we can set (Vi,Vi) '■= (l^n'Qn) wnere (£i,n,£i,n) are as m above 
proposition. Introducing a "macro" scale time step h = nho, the discrete 
model (|4.2p corresponds to the "aggregated" dynamics of (|4.9I) . Set for all 
(01,02) ip(0i,0 2 ) '-fpn (01, 02)- With the notations of Section O one 

derives that q n (Oi,02) = /n(0i;02) satisfies (A2) with the above ^. 

APPENDIX A: PROOF OF THE LOCAL LIMIT THEOREM O 

From now on, we use the following notations for multi-indices and powers. 
For v= (vi,...,v 2 d) G N M , = (xi, yi, y d )* set 

|z/| = Ui + ... + l/ 2c j, V\ = Vi\...V 2 d}-i 

(x,yy = x\\.. x^\.. tf*,ir = D»\...Dl* d D% + \..D2 d . 

A.l. Preliminary controls on the discrete kernel. We first give 
some controls for the kernel Hh(tj, (x,y), (x' ,y')). Namely, the following 
Lemma states that the difference between H^, '■= (L — L)ph and an 
additional remainder term is small, i.e. has the order announced in The- 
orem 14.11 

Lemma A.l (Control of the discrete kernel) For all j € set 

p 2 = tj . One has 

\{H h - K h - Mh)^,^),^' ,y'))\ 
(A.16) < Ch^p-'Ux' -x,y'-y- ^p 2 ). 

where Q p is as in Lemma \4-1\ and for j G (1,N\, 

K h (tj, (x,y), (x',y')) = (L - L)p h (tj, (x,y + xh), {x',y')), 

i.e. Kh is the difference of the generators associated to the initial and frozen 
diffusion processes between and tj applied to the density of the Markov 
chain with a slight change for the initial point in the y component, 

3 

(A.17) i\4fo,(z,y),(x',i/)) = £M£(t i ,(x,y),(x',y')), 

k=l 

where the (M^^g^^j are defined in Appendix^ 
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For j = 1 we set K h (t jt (x,y), {x',y')) = 0, 

M h (tj, (x,y), {x ,y')) = H h (tj, (x,y), (x',y')). 

The proof is postponed to Appendix [B] Prom this proof one also derives that 
the terms appearing in Lemma IA.1I are controlled with the following: 

Lemma A. 2 There exists a constant C s.t. for all < j < N , for all (x,y) 
and (x',y') in R 2d 

3 

(\K h \ + \M h \ +Y) M l\ + \H h \)(t j ,(x,y),(x', 1 /) 

i=l 

^ n -1a ( i i X + x' 2 \ 
<Cp Qp\x -x,y -y — p J , 

with C p as in Lemma \4-l\ Here again p = y/tj. 

The key fact is that the previous bound provides an integrable singularity in 
P- 

A. 2. Comparison of the parametrix expansions for the compen- 
sated diffusion and Markov chain. We first state an auxiliary result 
concerning the behavior of the iterated discrete kernel applied to the density 
of the frozen Markov chain. 

Lemma A. 3 There exists a constant C (that does not depend on (x,y) and 
(x',y')) such that, for all < j < N, re {0,j}, 

f I I X -\- X n\ 

x Xp I x - x, y -y — p I , 

where xp an d S' are as in Theorem \4-l\ 

To prove the lemma it is sufficient to repeat the proof of Lemma 3.11 in 
[KMOOj with obvious modifications concerning the additional arguments in 

y,y'- 

Lemma A. 4 For < j < N the following formula holds: 

p h (tj, (x, y), (x',y')) = Y.(p®h (M h + K h )^) (t jt (x, y), {x', y 1 )) + R, 

r=0 

where \R\ < Ch 1 / 2 p~ 1 Xp( x ' — x,y' — y - ^^-p 2 ) for some constant C. The 
function xp * s as in Theorem \4-l\ 



< 



L 



PARAMETRIX FOR SOME DEGENERATE PROCESSES 



23 



The proof follows from Lemmas 14 . 1 1 and [A . 2 1 and is analogous to the proof 
of Lemma 3.13. in |KM00j . □ 

Let us now compare the parametrix expansions of the compensated diffu- 
sion and Markov chain at time T. From Proposition 12.21 (|3 . 2[) and Stirling's 
asymptotic formula for the V function we have 

TV 

(A.18) p(T, (x, y), (x',y')) = £ (p ® ff«) (T, (x, y), (x>,y>)) + R u 

r=0 

where \Ri\ < Ch l /' 2 p c (T, x' - x, y' - y - ^T), with p c as in Theorem EH 
By Lemma lA.41 
(A.19) 



TV 

p h (T, (x, y), (x',y')) = £ (p ® h (M h + K h )W) (T, (x, y), (x',y')) + R 2 

r=0 

where 



\R 2 \ < Ch^T-^x^x' -x,y'-y- £±^T). 
Because of (|A.18I) and (|A.19I) . to prove the theorem it remains to show that 



IA 



TV 



TV TV \ 

E (p®H^)-J2(p®h (M h + K h )W)\(T, (x,y),(x',y')) 

r=0 r=0 J 



(A.20) < C(l + \x'\)h l l 2 XvT (x' -x,y'-y- —r—T). 



x + x' 



Note that |Ajv| < Si + 52 + S3 + ^4, where 

TV TV 



s 



So 



5, 



£ -J2 (p®hH^) (T, (x, y), (x , y')) 

\r=0 r=0 / 

TV TV \ 

r=0 r=0 / 

TV TV \ 

E(p% ffW ) -E(P®" (^ + #) (r) ) (r,(x,y),(x',y')) 

Vr=0 



V 



5, 



\r=Q 



E (p®* + - E (p®^ + ^) (r) ) ) ( T ' (sV)) 



V 



r=0 
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where H(t,(x,y),(x',y')) = H(t, (x,y + xh),(x' ,y')) is a shifted operator 
introduced for the comparison with K^, see the proof of Lemma [A. II in the 
AppendixlBlfor details. 
We shall show 

5 1 < C/.%( a /-i, ! /- ! ,-^T),iG{l,3 I 4} 1 

5 2 < C(l + |x , |)/ i 1 / 2 p c (^,( a; ,y),(x^y , ))■ 
This is done in Appendix ICl 

APPENDIX B: PROOF OF LEMMAS EH AND [Aj 

B.l. Proof of Lemma IA.ll For j = 1 we have p = y/h. By definition 
of H h 

H h (h, (x,y), {x ,y')) = p~ 2 (p h -p£ x ,y )(h, (x,y), (x',y')) 

Thus, recalling q n stands for the density of the random variables appearing 
in schemes (ED) , (l4~2l) 



\H h (h,(x,y),(x',y'))\=h-^ 

1 



1 



y/deta(x' , y' — x'h) 
where 

a~ 1 (x' , y' — x'h)(x' — x — b(x' ,y')h) 



y/deta(x, y) 
(u,v) 



(u + Si, v + 5 2 ) 



u 



a 1 (x,y)(x' - x - b(x,y)h) 
, u+di = 



-x, , 1 , u ,y' -y-{x + ±j n hb(x',y'))h 
v = a {x,y-xh) , 

, s -i f sV'-V- (x + ^ n hb(x,y))h 
v + 6 2 = a (x,y) 



h 3/2 



Note that 
B.21) 



1 



< C 



\J deta(x', y' — x'h) y/deta(x, y) 
(B.22) 

\q n (u + 5t, v + 5 2 ) - q n (u,v)\ < Ci/j (\u,u + 5i] 7 , [v,v + 5 2 ] 7 ) (|<5i| + \5 2 \) , 
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where [X, Y] := (1 — j)X + jY, 7 £ [0, 1], for X, Y matrices or vectors. One 
has 



[u,u + 8i] r 



rr 



~ x {x' ,y' - x'h),a 1 (x,y)] 7 (x / - x - b(x',y')h) 



[v,v + 62], 



+7cr~ 1 (x, y){b(x , y) - b(x, y))Vh 
= [a~ (x', y' — x'h), a^ 1 (x, y)] 7 0"(x', y' — x'h)u + jRi, 

l<T-\x>,y> - x>h),a-\x,y)] 7 (y> - y - {x + ^ x '/ )h )h) 

h 3 / 2 

+-ya~ 1 (x,y)(b(x' ,y') - b(x,y))^j-Vh 
:= [a~ 1 (x',y / - x'h), cr _1 (x, y)] 7 a(x', y' - x'h)v + 7-R2- 
Assumptions (B), (UE) also yield 

\Ri\ + \R2\<Ch^ 2 , 
3d,C 2 > 0, V(eR d , < 7 < 1, 

Ci|e| < \W~\x' ,y' - x'h)^- 1 ^^)]^' ,y' - x'h)ti\ < C 2 |£|. 
We also have 

\h\ + < chh (i + m 2 + m 2 ) . 

Recall that from (A2) and our choice for ip, < 1+ |( M ^^d+i ■ From 

(|B.21|) . (IB.22|) and the above computations we get 



\H h (h,(x,y),(x',y'))\ < Ch~ l > 2 h~ 



2d 



1 + Id 2 + \v\ 2 



(l + \(u,v)\ s + 2d + 1 
(B.23) < Cp^C^x' -x,y> -y 



■~y» I rv* f 



For 1 < j < N, we proceed like in the proof of Lemma 3.9 in [KMOOj . We 
get that 

Hh(tj, (x,y), (x',y')) = (H h - H h )(tj,(x,y), {x ,y')) 

where 

H h (t j ,(x,y),(x',y')) = h- 1 f q n {9 1: 9 2 ) x 
(B.24) \\(x + 7 X (0i), y + xh + rf (0 2 )) - A(x, y + xh)\ dBM, 
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H h (t j ,(x,y),(x',y'))=h- 1 [ 9 „(0i,0 2 ) * 
(B.25) [A(x + 'y 1 {9\),y + x/i + 7 2 (#2)) - A(s, y + x/i)] d0id0 2 , 

with A(u, v) = p h (tj - h, (u,v), (x',y')), 

f{9 x ) = hb(x, y) + Vha(x, y)9x, l\9 2 ) = j^M^ + y)fl 2 j fc, 



and 



j 1 (9 1 ) = hb(x',y') + Vha(x',y' -x'p 2 )9 u 



+ Vha{x\ y' — x'p 2 )9 2 ^ 



Using a Taylor expansion at order three for A in (|B.24p and (IB.25j) we 
obtain 

H h (tj, (x,y), (x',y')) = j^Tr(a(x,y) - a(x',y' - x' p 2 ))D 2 x \{x, y + xh) 
+(&(x, y) - 6(x', y'), V x A(x, y + xh) + -^-X7 y \(x, y + x/i)) 

+|^ ((D 2 x X(x, y + xfc)&(a:, y), b(x, y)) - (i^A(x, y + x/i)6(x', y'), v'))) 

+ y (tr(^A(x, y + xh)(a(x, y) - a(x', y' - p 2 x')))) ^ 7n (l + 

+^f- ((D 2 y \{x, y + xh)b(x, y), b(x, y))- (D 2 y X(x, y + xh)b(x',y'), b(x', y')) 



{h- 1 I d9 1 d9 2 dq n (9 1 ,9 2 ) 
(Dl x \(x,y + xh)j 1 (9 l ),j 2 (e 2 )) - (Dl x X(x,y + xh)f(9 1 ),y 2 (9 2 )))} 



+3/T 1 J de i d0 2 £ dS(l-S) 2 q n (9 1 ,9 2 



\u\=3 



D v X{x + Sf{e x ),y + xh + 5j 2 (6 2 )) 

(i 1 (e 1 )^ 2 {9 2 )y 



3h^ J Mid9 2 J dS(l-S) 2 q n (9 1 ,9 2 



kl=3 



(B.26) 



xD"\(x + ^(^i), y + xh + <57 2 (0 2 )) 
:=/ + // + m + IV -V, 
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where we denote D 2 X(x,y + xh) (resp. D 2 X(x,y + xh), D 2 x X(x,y + xh)) 
the R d (g> M d matrices (d x . jXj X{x,y + xh)) ^^^2 (resp. {d yi>yj X{x,y + 
x/l ))(i,j)6[M] 2 ' (d Vi , Xj \(x,y + x/i)) (ij -) g [ 1>d ]2). 

In the sequel, a useful result is the following. There exists C > s.t. for 
multi-indices a,/?, |a| < 3, \/3\ < 3, 

|d^A(x,y + xh)| < Cp-^l+^^DCp^'-x^'-y 

(B.27) < Cp-^+^CpjV-^y'-y 

This assertion can be proved similarly to Lemma 3.7 in |KM00j . 
Note now that 

I = (L — L)p h (tj, (x, y + xh), (x',y')) 

+\~Y-{K x i y) - K x 'i y')^yK x i v + xh )) 

+{L - L)(X(x, y + xh)- p h (tj, (x, y + xh), {x',y'))) j 
:= {K h + Ml)(t 3 ,{x,y),{x',y')). 
From the above equation and (|B.27|) we get 

\Mk(tj, (x,y), (x',y'))\ < Cp^ 1 Q p {x - x, y - y - ^—^-p 2 ). 

(B.28) 

Using similarly (|B.27|) and tedious but elementary calculations, one can 
split in II, III the terms that give the expected order, i.e. bounded by 
C\fhp~ x Q p (x' -x,y' -y-^-p 2 ) and denoted below by R\{tj, (x,y), (x',y')), 
and those that give an integrable singularity in time, i.e. bounded by Cp~ 1 C p (x' — 
x,y' — y— ^-p 2 ) and denoted below by M%(tj, (x,y), (x',y')). 



— xh 



x + x' 
2 
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It remains to estimate IV — V in (|B.26l) . To this end write, 

iv-v = 3h~ 1 ^Jde l dd 2 j\5{i-5) 2 q n {e 1 ,e 2 )l 

|i/|=3 

{{i 1 {e 1 )^ 2 {e 2 )Y - (f{e l )^ 2 {9 2 )Y)D"\{x + af^), y + xh + sf(9 2 )) 

+(7 1 (#i), f(0 2 , ))" E tdaD^\{x + Sfi^) + a^ 1 - f)(0i), 

y + x/i + 57 2 (^ 2 ) + a<5(7 2 (# 2 ) - 7 2 (^))) 
(5(7 X - 7 1 )(^i), 6tf(9 2 ) - f(0 2 ))Y }:= M^tj, (x, y), {x\y')). 

Computations involving (|B.27|) yield 

(B.29) \Ml{tj, (x, y), (x>, y'))\ < Cp-\ p {x' -x,y'-y- ^- P 2 ). 

We refer to the proof of (3.80) p. 584 in |KM00] and Appendix El for 
additional details. This completes the proof. □ 
The proof of Lemma IA.2I then follows from the previous proof, (lB~27l) . 
(|B~28]) . (lB~29l) and (lB~26l) for j £ (1, Nj and (lB~23l) for j = 1. 

APPENDIX C: CONTROL OF THE (S/)/ e [i, 4 ] 
C.l. Control of Si. Set 

oo 

Pd(T, (x,y), (x',y')) = ^p® h # (r) (T, (x,y), (x',y')). 

r=0 

From Proposition 12.21 one has 

{p-p d )(T,(x,y),{x',y')) = (p <g> H - p ® h H)(T, (x, y), (x , y')) 

+{P-Vd) ®h H(T, (x,?/), Obl- 
iterating the previous identity we get 

(P - Pd){T, (x, y), {x',y')) = (p g> H - p ® h H) ® h <p(T, (x, y), (x',y')), 
(C.30) 

where Vj £ [0, JV — 1], V(u,t>) £ M M , 

oo 

¥>(T - (u, u), (s', y')) = X>£ r) (T " *i> ( u ' «). (^^))- 
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Let us first give a bound for Pj(u, v) := (p<g)ii"— p®hH)(tj^ (x,y), (u, v)), j £ 
[0, iV], (it, v) G IR M . First, from the previous definitions of the continuous 
and discrete convolution operators, Pq(u,v) = 0, in the sense of generalized 
functions. For j > 1 write 



3-1 



Pj(u,v) = [ dt [ dwdz\ u>v ) (t, (w, z)) - A( U) „) (t», (w, z)), 

\u,v)(t,( w , z )) P(t,(x,y),(w,z))H(tj -t,(w,z),(u,v)). 
A first order Taylor expansion and Fubini's theorem give 

rU+i /•! . 



Pj(u,v) 


E/" 


Qi(u,v,s) 


:= / i 


rpO 

3 


:= / dt 
Jo 



(C.31) x (H(tj — t, (w, z), (u, v)) — H(tj, (x, y), (u, v))). 

From Lemma 13.11 Theorem 12.11 we obtain 

Tf < CVhp c (tj,(x,y),(u,v)), 
\d s X {u>v) (s,w,z)\ < Cis-^tj-sy^ + itj-s)- 3 / 2 ) 

xp c (s, (x,y), (w,z))p c {tj - s, (w,z), (u,v)). 

Plug now the above control in (|C.31j) . we get 

/ L0-i)/2j 

Pj(u,v) < Cp c (t j ,( x , y ),(u,v))(h l / 2 + h 2 [tj 1/2 £ ti 1 



1=1 



+t] 1 E (^-^ + i)- 1/2 + eV 3/2 

i=L(i-l)/2j+l i=l 

< Ch^ 2 p c {tj,(x,y),(u,v)). 



Hence, from (IC.3Q|) and a suitable version of (|3.2p for the discrete convolution 
operator we derive 

\{p-p d )(T,(x,y),(x',y'))\ < Ch^p^T, (x,y), (x' ,y')). 



The bound for Si can be derived using once again (13.21) for both the contin- 
uous and discrete convolution operators and the asymptotics of the Gamma 
function. 
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_ C.2. Control of S 2 . Define for r G [0,iV], T r := (p ® #M -p® 
H^)(T, (x,y), (x' ,y')). For r = 1, with the notations of Lemma l3.1l one gets 

JV-l f i . 

\Ti\ < Ch 2 J2 d5 /tJ^CT-^r^HPcCti.Ca:.!/), («,«)) 

xp c (T — tj, (u, v + (x ,y'))dudv. 

Also, for a different constant c/ than the one appearing in p c , V? G [1, iV — 1], 

p c (T - t jt (u, v + 8uh), (x', y')) < C(T - tj r (3k+d)/2 

( ,\\x'-u\ 2 W -v-*¥(T + 8h-U)\ 2 \\ 

(C.32) < Cp c ,(T + 5h-t v (u,v),(x',y')). 

Hence, up to another suitable modification of the constant in order to have 
the semigroup property 

N-l i 

|Ti| < Ch 2y £tJ 3/2 (l + (T-t j )- 1 ^\a/\) d5p c (T + 5h,(x,y),(x',y>)) 

< Ch l ' 2 {\ + T^Vl) f 1 dSp c (T + Sh, (x, y), (x', y')). 

Jo 

Write now, for all r > 2, 

T r = p® h H( r ~V ® h (H - H)(T,(x,y),(x',y')) 

+(p® h -p® h H^) ® h H(T, (x,y), (x',y')) := T rl +T r2 . 

The term T r \ can be handled as T\ exploiting the control 

\d z p® h H^(t,(x,y),(w,z))\ < C r+1 t^/ 2 p c (t,(x,y),(w,z)) 

xp((l + i)/2,l/2). 

i=0 

For T r 2 one uses the control of step (r — 1). Completing the induction one 
derives 

|S 2 |<Cfc 1 / 2 (l + |s'|) sup p c (T(l + 6),(x,y),(x',y')). 

56 [0,1] 

Note that this term is the only one for which we have a linear contribution of 
the terminal variable. This is, because of the shift, in some sense unavoidable. 
Also, the previous trick in (1C.32|) adds the constraint to take a supremum 
w.r.t. to a twice larger time interval as the initial one. 



PARAMETRIX FOR SOME DEGENERATE PROCESSES 



31 



C.3. Control of S3. For r = 1 we have to control 

p® h M h (T, (x,y),(x',y')) = 

3 N-l 

J2 h Yl / dudvp(tj, (x,y), (u,v))M l h (T - tj, (u,v), (x',y')) 



3 N-2 
i=l j=0 

The term Wjv-i needs to be handled by a different technique than the 
other ones. Write 



Mn-i = J dudvp(T — h, (x,y), (u,v)) 
x (Ph-Ph) (h, (u,v), {x',y')). 

Set V = ( (r"fa)i/2 1 V (T-h)^ 2 Write now \u — x\ = \x' — x + u — x'\ 



v - y — (T - h) 



ry I ry rf r*J 

y — y — T H — h + v — y + uh 



I ,_ y> y x + x T I X ~ X k \ 

Set U = 77fr-rw72 1 (T-h\3/2 2 — ) V := U+R. Recall also from Lemma 



(T-/l)V2 > (T-/l)3/ 2 

Othatp < Cp c . Hence, for all Z G N*, 3C := C(Z), p(T-h, (x,y), (u,v)) < 
(T - h)~ 2d l+ yv\z ■ From the basic identity 1+ ^ +R ^ < max ( 2 i ^+(^l R D ) an d 
the definitions of the models (14. ip and (|4.2|> . using the same techniques as 
in the proof of Lemma fA. II for the case j = 1 one gets: 

< 1 + m z j du'dv'(l + \( U ',v')\ z M-u', -v') 

< Ch^C^ix'-x^'-y-^T), 



taking Z = S — 4 for the last inequality. 

Also, from the definitions of the (M^)j g j 1)3 j in the previous section and 
using freely its notations, we derive for all j G [0, N — 2]: 

\Ml{T-ti,{u,vW^))\ < h(T-t j r^ 2 C p (x'- U ,y'-v- { -^±^(T- 
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N-2 x + X 1 

from which one gets h ^ \h,j\ < Ch 1 ^ 2 C i ^?p(x' — x,y' — y — T). The 

j=o 2 
terms in M% coming from II in (lT3~26l) can be handled as M\. For those 
coming from III, i.e. crossed derivatives, the contribution associated to j = 
is easily analyzed and for j > 1 an integration by part w.r.t. u leads to the 
same control. The trickiest term to analyze is M^. Exploiting thoroughly 
(|B,27|) and Lemma I3~T1 the proof is similar to the one in |KM00| . see p. 578 
control of (3.45), that relies on suitable integration by parts. We omit the 
details here. Actually, for r > 1 it can be shown by induction that 

(p ® h -p® h (H + M)«) (T, (x, y), (x 1 , y')) 
h l / 2 C r+1 , . , x + x' s 

- Twrm XVT(x - x ^- y - — n 

which gives the control. 

C.4. Control of 54. One can show that Lemma 14.11 is still valid for 
the derivatives of the frozen densities. Using this result and Lemma [A. 21 the 
proof is then similar to the one of |KM00| . 
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